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Abstract 

Precise magnetic structures of RMn20t,, with R = Y, Ho, Bi in the commensurate/ferroelectric 
regime, have been determined by single-crystal neutron diffraction. For each system, the integrated 
intensities of a large number of independent magnetic Bragg reflections have been measured, allow- 
ing unconstrained least-squares refinement of the structures. The analysis confirms the previously 
reported magnetic configuration in the ab-plane, in particular the existence of zig-zag antiferro- 
magnetic chains. For the Y and Ho compounds additional weak magnetic components parallel to 
the c-axis were detected which are modulated in phase quadrature with the a-b components. This 
component is extremely small in the BiMn205 sample, therefore supporting symmetric exchange 
as the principal mechanism inducing ferroelectricity. For HoMn205, a magnetic ordering of the Ho 
moments was observed, which is consistent with a super-exchange interaction through the oxygens. 
For all three compounds, the point symmetry in the magnetically ordered state is m2m, allowing 
the polar 6-axis found experimentally. 



PACS numbers: 25.40.Dn, 75.25.+Z, 77.80.-e 



I. INTRODUCTION 



Materials in which there is strong interplay between ferroelectric and magnetic ordering 
are attracting current interest due to the possibility of controlling the electric polarization 
by application of a magnetic field or conversely the magnetisation with an electric field. 
In the last few years, studies have been focussed on systems, such as TbMnOs and 
i?Mn205[3], for which the ferroelectric transition coincides with the transition to a com- 
plex antiferromagnetic ordered phase at low-temperature. These systems are classified as 
improper ferroelectrics, since the primary order parameter is the magnetization, and their 
electrical polarization is at least an order of magnitude lower than in proper ferroelectrics ^ . 
Nonetheless, the coupling between the magnetization and the electrical polarization is strong, 
providing an exciting playground in which to study the microscopic mechanisms that govern 
the magneto-electric interaction. Moreover, the aforementioned systems commonly present 
a high degree of magnetic frustration, which seems closely linked to the appearance of 
ferroelectricity [sl . 

The nature of the low-temperature state in these materials has been investigated by a va- 
riety of experimental techniques sensitive to small ionic displacements in a crystal (Raman, 
IR spectroscopy) and scatterin g te chniques sensitive to the magnetic order (Neutron, 
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12l |. An increasing number of theoretical models 



resonant X-ray scattering 
have also appeared over the last two years that discuss the origin of the magneto-electric 
coupling!^, 1^, Q, 15, IgI. Two microscopic mechanisms have been put forward: on one 
hand the antisymmetric inverse Dzyaloshinskii-Moriya (DM) interaction, requiring a non- 
coUinear magnetic arrangement and on the other hand exchange-striction due to symmetric 
exchange. In the DM model, the polarization direction P (atomic displacement direction) 
is given by the expression: 

P oc ei2 X (Si X ^2) (1) 

where Si and S2 are neighboring spins and ei2 is the propagation direction. The propor- 
tionality constant is related to the strength of the spin-orbit coupling. It seems that this 
mechanism explains the origin of ferroelectricity in many systems recently investigated such 



as TbMnO^ 
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15| . Exchange-striction, which is intrinsically stronger that DM, does not 



require non-coUinearity, and has been proposed to explain the appearance of ferroelectricity 
in RMuoOpiflGl and in the 'E' -phase of HoMnOa 17 1. Recent work by Infra- Red absorption 
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confirmed that both mechanisms described above are relevant, as suggested by different po- 
larization selection rules for the electromagnons in RMnO^ and /?Mn2 05[7]. It is important 
to note, that due to the extremely small shifts of atoms in the ferroelectric state with respect 
to their positions in the centro-symmetric paraelectric phase, these displacement patterns 
have not yet been properly determined. The mechanism which drives the ferroelectric tran- 
sition can therefore often be more easily inferred from analysis of the magnetic configuration 
derived from neutron diffraction experiments. 

There has been continuing interest in the complex magnetic structures of i?Mn2 05 with 
R = Rare earth (space group Pbam) since the 1960s Q, Q, Q- 

These insulators order 

antiferromagnetically at low temperature (T^v ~ 40 K) with a propagation vector (k^, 0, k^) 
and on further cooling undergo a series of magnetic transitions to both commensurate and 
incommensurate magnetic phases. These materials become ferroelectric at a temperature 
slightly below the Neel transition at T^v with values of electric polarization, ranging from 
20 to 100 nC.cm"^ depending on the R ion, being largest in the commensurate phase. The 
component of the magnetic propagation vector along c (k^) depends on the size of the R 
cation, which suggests that the magnetic exchange along the c-direction is extremely sen- 
sitive to the interlayer couplingjo], U^- Bi with r/=131 pm has the largest ionic radius of 
the series and in BiMuaOs k^=l/2 whereas k^=l/4 for YMuaOs, HoMusOs and TbMuaOs 



with r/=116,116,118 pm respectively 



32] ■ On the other hand, the onset and strength of 



ferroelectricity does not seem to be dependent on tj. 

The magnetic structures of this class of materials have mostly been derived from neutron 
diffraction on polycrystalline samples. Recently however, Noda and co-workers have re- 
ported results obtained by single-crystal neutron diffraction for the commensurate magnetic 
phase of YMn205j8]. Their analysis suggested the presence of a small c-axis component 
of magnetic moment on the Mn sites, which had not been detected in previous powder 
measurements jl^. This additional component is in phase quadrature with the major a-b 
component and introduces a cycloidal modulation of the magnetic moments which could, in 
principle, induce ferroelectric order through the DM interaction. 

In the light of this recent work, we have undertaken determination of the commensurate 
magnetic structures of three compounds BiMuaOs, YMuaOs and HoMuaOs using single- 
crystal neutron diffraction. Variation of the cations allows a comparative study in which 
their influence on various characteristics of the magnetic structures can be distinguished. 



3 



The aim of the study is twofold: first to search systematically for the presence of the c-axis 
modulation in several members of the series so as to gauge its importance in promoting 
ferroelectricity. Determination of the magnitude of this component for all three compounds 
is especially important because the magnetic propagation vector in BiMn205 (k^=|) dif- 
fers from that in HoMn205 and YMn205 (kz=^). Secondly, by studying HoMn205, we can 
investigate the influence of a magnetic R ion on the arrangement of the Mn moments and 
determine the magnetic ordering of the Ho sublattice itself. Comparison with YMn205, 
which has the same propagation vector but a non- magnetic R site is of particular interest, 
since the values of the electric polarization in the commensurate regime of the two systems 
are significantly different. 

II. EXPERIMENT 

Single crystals of i?Mn205 {R—Bi, Y, Ho) were grown using B203/PbO/PbF2 flux in a 
Pt crucible. The flux was held at 1,280°C for 15 hours and slowly cooled down to 950°C 
at a rate of 1°C per hour. Crystals grew in the form of cubes. The samples used for 
the present work were respectively of sizes: ~ 4x4x4 mm^ for YMn205, 2x2x2 mm^ for 
HoMn205 and 1x1x3 mm^ for BiMn205, where the three dimensions refer to the [110], 
[TlO] and [001] directions. For the HoMn205 crystal, additional faces (100) and (010) were 
also visible. Single crystal diffraction measurements were performed on the four-circle 
diffractometer DIO at the ILL (Grenoble, Prance). Samples were checked for quality and 
pre-aligned with the (001) direction oriented along the vertical axis using the OrientExpress 
facility (ILL). Samples were mounted on standard aluminum pins in a He exchange gas 
cryostat. Data were collected with an incident neutron wavelength of A = 2.36 A using an 
80 mm^ 2D micro-strip detector. Q scans around each Bragg reflection were performed, 
with variable counting time depending on the intensity. Peak integration was performed 
using the program Racer (ILL), in two steps. First, a library was built by fitting ellipsoidal 
shapes to a set of strong reflections (I>3(t), these shapes were used in a second pass to 
integrate all reflections. Data were collected for all three samples in the ordered magnetic 
commensurate phases. An additional datasct was collected for the YMn205 sample in the 
low-temperature incommensurate phase, but the results will be reported elsewhere. The 
data sets for YMn205 were collected at T = 25 K (nuclear and magnetic scattering), those 
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for BiMn205 at T = 10 K (nuclear and magnetic scattering) and for HoMn205 at T = 
50 K (nuclear scattering) and T = 25 K (magnetic scattering). For each data-set, the 
list of integrated intensities obtained were corrected for Lorentz-factors and normalized to 
an arbitrary monitor count. Determination of both the nuclear and magnetic structures 
were carried out using the FullProf program suite • For each sample the following 
procedure was applied: the published crystallographic models, space group Pbam were 
used to refine the nuclear data-sets. Because a limited Q-range is accessible on DIO, the 
crystallographic parameters were fixed, with the exception of a global thermal parameter to 
account for different temperatures. In addition, a global scale factor was refined together 
with six parameters for extinction correction, following the formulation of Becker-Coppens 
for anisotropic extinction (extinction model number 4 in FullProf). This correction is 
particularly important for YMn205 because of the large size of the crystal. The scale 
factors obtained by refining the nuclear structure were used and fixed for refinement of the 
magnetic structures. 



III. RESULTS 

A. Crystal structure and formalism 

The crystal 



detail in 
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structure and topology of the magnetic interactions have been described in 
23(1 . The different atomic positions for Mn^+, Mn^"*" (and Ho^"'') have been 
numbered according to reference {23 1. To facilitate the identification of all magnetic sites 
in the figures, additional labels have been included: a for Mn^"*", b for Mn''"'"(c for Ho^+). 
The crystal structure of i?Mn205 (isostructural throughout the series) is displayed in Fig. 
[2] and the correspondence between labels and atomic positions can be found in table HTl 
(a) is five-coordinated by oxygen ions, in square-pyramid geometry, while Mn^"*" (b) 
is six-coordinated in distorted octahedral geometry. There are four symmetry equivalent 
Mn^+ sites and four Mn^+ sites per unit-cell. The pyramids centered at a2 and a3 (also 
al and a4) share an edge of their basal plane. All Mn^+ are located on a mirror plane at 
z=0.5. The octahedra coordinating the Mn^+ sites bl and b2 (also b3 and b4) share edges 
to form chains running along the c-direction. The bl(b3) atoms with z~0.25 and b2(b4) 
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with(l-z)~0.75 lie below and above the Mn^+ layer. It should be noted that the sites bi and 
b2 also ba and are related by a center of symmetry. The Mn'^"'" octahedra share corners 
with Mn^+Os pyramids to form layers in the a6-plane. The R ions (labeled c for HoMn205) 
occupy the (x,y,0) site, of multiplicity four, and are all positioned in the z=0 layer separating 
two Mn^+/ Mn^'^ layers. 

The magnetic structures have been analyzed using the propagation-vector formalism, in 
which magnetic moments for all sites described above are expanded in Fourier series. We 
briefly describe the conventions used in the next paragraph. For a magnetic structure with 
a single propagation vector k , which is the relevant case in this study, a magnetic moment 
of atom type j in the crystallographic unit-cell 1 (m/j) is written: 



S^j ■ e-2--(^-^+'^^) (2) 



where the sum runs over k and - k if these vectors are not related by a reciprocal lattice 
vector (case of YMn205 and HoMn205) or uniquely k otherwise (case of BiMn205). Ri is a 
pure lattice translation and (pj is a phase factor. The magnetic structure is fully determined 
by giving a set of Fourier coefficients and phases {Skj,4>j) for each j, i.e seven components 
per site in the most general case since the Fourier coefficients are complex quantities. We 



note that because rriij is a real quantity, the constraint S-kj=Slj (where * denotes the 
complex conjugate) is required. It is important to mention that (pj are relative phases and 
are equal for sites belonging to the same orbit, i.e. sites that can be transformed into each 
other by an operation of the magnetic little group. It is different than a global phase, 
that can't be determined by diffraction, which rephase the entire magnetic structure, i.e. 
apply simultaneously to all magnetic sites. Additional constraints exist between {Skj,4'j) of 
symmetry-related sites. These constraints are determined by representation analysis and 



will be briefly described in the following sections, based on previously published work 24]. 
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B. YMnaOs 



1. Magnetic structure 
The propagation vector for t 



uhe Y-compound in the commensurate phase (CM) is =(1/2, 
0, 1/4). Symmetry analysis 2J] indicates that the Mn^"*" sites is spht into two orbits because 
the rriz mirror-plane operation does not belong to the little group. The first and second or- 
bits contain respectively the bl/b3 sites and b2/b4 sites. The (pj phases for each orbit can 
therefore be non-equal and treated as free parameters during the refinement. On the other 
hand, the phases of the Mn^+ sites, belonging to a single orbit, must be equal. We arbitrarily 
chose to fix the phase of the Mn^"*" sites at 0.125 (one phase must be fixed since structure 
factors do not depend on a global phase), for comparison with results reported earlier from 
powder data [l^ describing a configuration with equal moments. In this case, representa- 
tion analysis does not impose any constraints on the magnitude and direction of the Fourier 
coefficients for both Mn^+ and Mn^^ sites|24!]. For each site, six independent Fourier coef- 
ficients are in principle refinable. The unconstrained model has therefore fifty parameters 
and describe the most general magnetic structure. For typical magnetic structures described 
by either spin-density waves (SDW) or simple cycloidal modulations, some of these param- 
eters will be fully correlated. Refining an unconstrained model against the 355 independent 
magnetic reflections fails due to the large number of parameters and the presence of strong 
correlations. In order to find possible models, a global optimization algorithm, using the 



simulated annealing procedure described in 



25| . was employed. Systematically, the Fourier 



coefficients for all magnetic sites were found to be principally along the a-axis, whether real 
or imaginary. In addition, a weaker 6-axis component was found with the same character 
(real or imaginary) as the a-component. These two components define a vector, say vi in 
the aft-plane. A small c-axis component in phase quadrature with the a and 6-components 
was also identified. The structure therefore has a small helicoidal modulation, with rotation 
axis given by the cross product of vi with a vector along z. A model with this configuration, 
i.e one in which the Fourier coefficients for all magnetic sites are real in the aft-plane and 
imaginary along c, has been constructed, it has 26 refinable parameters. The lower number 
of free parameters allowed a direct least-square refinement to be made which converged after 
a few cycles. The structure factors calculated with this model are shown, plotted against 
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those observed, in Fig. [T] The model was found to reproduce the observed data very well 
with a Magnetic R(F) factor [33| of 4.57% (Fig. [1]). For this particular sample, the large 
value of x^=163 is explained by high statistics of the data obtained on a large sample with 
long counting time. The refined values of the Fourier coefficients are displayed in Table [TTl 
and the magnetic structure projected along two different crystallographic directions in Fig. 
[3] and Fig. |H The magnetic arrangement in the a6-plane (Fig. [3]) is equivalent to that 



derived in our previous work from powder data[16|. The moment directions are within ~10° 
of the a-axis for Mn^+and ~14° for Mn'*+, defining zig-zag antiferromagnetic (AFM) chains 
running along a. One of these chains links bl, b2, a2 and a3, the other b3, b4, a4 and al. 
In addition, there is a significant c-component on all magnetic sites as suggested by Noda et 
al.j^, which was not found in our earlier analysis of powder data. The c-components have 
almost the same magnitude as the 6-components on both Mn sites. Their imaginary char- 
acter means that the modulation is in quadrature with the real in-plane components. This 
can be seen in Fig. |H for example for the a2 site the in-plane components of the magnetic 
moment follow a (- H — h -) sequence across four unit-cells along c whilst the c-component 

follows a ( h +) sequence. This result is also in agreement with the work of Noda et al. 

The refined values of the phases for sites bl/b3 and b2/b4, respectively 0.095(2) and 0.156(2), 
are slightly different from the fixed values of 0.125 used for fitting powder data. Conse- 
quently, and in addition to the existence of a c- component, the magnetic moments on site 
bl and b2 (as well as b3 and b4) are of different magnitudes as shown in FigHl We also note 
that there is a symmetric deviation of these phases with respect to Mn^+ phase of 0.125 ( 
(0.125-5) for the first orbit and (0.125+5) for the second orbit). The identification of a small 
c-component has direct implications on the nature of the magnetic configuration, previously 
described as a pure SDW: small cycloidal modulations propagating along the edge-sharing 
Mn^+Oe chains are observed in the ac and he -planes (Fig. [5]). In the ac-plane, the moment 
rotation direction is clockwise for Mn^+ chains formed by the bl and b2 sites (left panel of 
Fig. HI) and anti-clockwise for chains formed by the b3 and b4 sites (right panel of Fig. Hj). 
In the 6c-plane, the weak cycloidal modulation, highlighted in Fig. [Sj rotates in the same 
direction along bl/b2 chains and b3/b4 chains. The Mn^+ moments are also slightly tilted 
along the c-direction, remaining almost collinear to adjacent Mn'^+ magnetic moments in 
a given AFM chain. Therefore, the presence of antiferromagnetic chains running along the 
a-axis, identified earlier as an important characteristic of these systems, remain valid despite 
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the small c-component. 



2. Domains 

The arm of the star of k is made only of the two vectors ^ = (|,0,|) and - k since the 
vectors (-|,0,|) and (|,0,-|) are related to k and -k by a reciprocal lattice vector. Therefore, 
a single k-domain contributes to the scattering. On the other hand, one must take into 
account the "orientation" domains, i.e. domains that are obtained by applying symmetry 
operations of the paramagnetic group that are not valid operations (lost) of the magnetic 
structure. Four symmetry operations of the paramagnetic group (combined or not with 
complex conjugation) leaves the magnetic structure unchanged, as shown in Table [B We 
refer to this domain as domain one (Dl). The other four symmetry operations all transform 
the magnetic structure described in the previous section into an inequivalent one (domain 
2, D2, shown in Fig. [6]). No other domains exist, and Dl and D2 are simply related by 
inversion symmetry. Inversion symmetry operates differently on different magnetic sites: for 
ba, b4, ai and a^, belonging to one of the zig-zag AFM chain, the application of inversion 
symmetry leaves the a and b moment components unchanged while it changes the sign of 
the c-component. In contrast for sites bi, b2, a2 and (the other AFM chain), inversion 
symmetry reverses the sign of the a and b components while preserving the c component. 
Since Dl and D2 are related by inversion symmetry, the diffracted intensities arising from 
both domains are identical for a non-polarized diffraction experiment and refining the data 
assuming any domain leads to the same result. We note that the two domains are also 
related by a simple rotation by 180° along the z-axis or the a;-axis. 

From the list of symmetry operations (Table. [T]) that leave the magnetic structure in- 
variant (modulo complex conjugation operation), one derives the point group m2m in the 
magnetically ordered phase, in agreement with that obtained earlier from co-representation 
analysis {2^. This point group is compatible with ferroelectricity along the fe-axis, irrespec- 
tive of the microscopic magneto-electric mechanism. 
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C. HoMnsOs 



HoMn205 orders AFM at T = 44 K first with an incommensurate propagation vector 
fc=(0.48, 0, 0.245), followed by a lock-in at T = 38 K to a commensurate =(|, 0, |). 
The magnetic structure has been determined in the commensurate phase from 381 inde- 
pendent magnetic reflections, using the same procedure as that described in the case of 
YMn205. However in this case, the magnetic ordering of Ho is an additional complexity 
which increases the number of free parameters. Attempting to refine the data without a 
magnetic contribution on the Ho sites, leads to poor agreement factors (R(F)=20.2%) as 
already reported from analyzing powder diffraction patterns 23|]. Preliminary measurements 
with resonant X-ray scattering at the Lm edge directly confirm the magnetic ordering of 
Ho. Symmetry analysis indicates that the four equivalent Ho positions belong to the same 
orbit and therefore only one phase parameter is required. However, there are no constraints 
between components of the moment of different sites, which were refined independently. 
The final refinement contained 45 parameters and was found to reproduce the data very 
well with structure factor R(F)=3.96% and x^=5.82. The magnetic structure is similar to 
the YMn205 one, with the appearance of almost collinear ...-Mn^"'"-Mn^"'"-Mn^+-Mn^+-... 
AFM zig-zag chains in the a6-plane. These chains make an angle of ~14° with respect to 
the a-axis. We again found a small c-component (Fig. [8]) of moment on all Mn magnetic 
sites, in quadrature with the in-plane components and of magnitude similar to that observed 
for YMn205. The refined phase for the Ho sites is shifted by 0.126(2) with respect to the 
Mn^+ sites. The Ho moments at the c3/c4 sites are parallel to the a-axis and those at the 
cl/c2 sites to the 6-axis (Fig. [7] and Table. IIIII) . It is possible that a smaller, tt/2 out of 
phase component along c may exist, but the refined magnitudes are only three times larger 
than the error bars. Although the magnetic configuration of the Ho sites is identical to 
that derived earlier by simulated annealing using powder data23|], the orientation of the 
Ho moments in the a6-plane is quite different. The values of the Ho moments, ranging 
from 1.1 to 1.3 fiB, are much smaller than the fully saturated value for a J=8 state. This 
can be partly explained by splitting of the fully degenerate J=8 state by the crystal field, 
but typical crystal-field energies in oxides will lead to a moment much larger than observed 



here 



261]. Additional reduction can be due thermal excitation, at 25K, of levels above the 



ground state (Ho^^ is a non-Kramer ion characterized by singlet levels). Moreover, dipolar 
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field calculations show no direct correlation between the local field generated by neighbors 
magnetic Mn sites and the magnitude and direction of the Ho moments. Therefore we argue 
that the Ho moment of ^1.1^ b in HoMn205 is mainly due to super-exchange interactions 
via the 2p oxygen orbitals. Exactly half of the Ho layers show ordered magnetic moments 
in the a6-plane, while the other half show only a small moment along the c-direction (Fig. 
|8]). Accordingly, in the former type of layers, we can clearly show that the directions of the 
moment are systematically arranged along the resultant of AFM super-exchange interactions 
with the four nearest-neighbor Mn^"*" sites (Fig. [7] and Fig. [8]). The net moment on each Ho 
sites is pointing either along the a-axis (sites c3 and c4) or 6-axis (sites cl and c2), resulting 
from the cancelation of respectively the b- and a-components of neighbor Mn^+ moments. 
Also, the resultant of the super-exchange is larger for sites c3 and c4 (along the a-direction) 
than for sites cl and c2 (along b), which should lead to different Ho moments on these 
sites, as we appear to observe experimentally (Table [In|) . A more precise determination of 
the respective magnitude of the moments on sites cl/c2 and c3/c4 is currently underway 
by resonant X-ray scattering. Calculations show that azimuthal scans on magnetic Bragg 
peaks at the Ho Lm edge are extremely sensitive to these parameters (as well as being Ho 
selective). 

The magnetic configuration of the Ho sublattice also breaks inversion symmetry since 
moments on sites c2 and c3 (respectively cl and c4) are pointing in different directions. As 
in the case of YMn205, there are only two "orientation" domains to consider, related by 
inversion symmetry and therefore contributing equally to the diffracted intensities. The sym- 
metry operations and corresponding domains listed in Table [T] are also valid for HoMn205. 
The magnetic structure of domain 2 including the transformed Ho moments (not shown) is 
easily derived from domain 1 by applying inversion symmetry. In this case, the magnetic 
point group is also m2m, which supports ferroelectricity along the fe-axis. 

D. BiMn205 

BiMn205 is the only member of the series that shows magnetic ordering below T=39 K 
with a propagation vector k =(|, 0, |) and no transitions to incommensurate magnetic order. 
Because of the commensurate propagation vector representation analysis leads to stricter 
constraints than those applicable to YMn205 and HoMn205, they are described in detail by 
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Munoz et al.|2l!]. No phases are required at any of the sites since k is half a reciprocal lattice 
vector. The magnetic representation contains two irreducible representations of dimension 
2. Only magnetic modes belonging to the representation Fi fit our single crystal data, 
in agreement with the results of Munoz et al. The magnetic moments on Mn^"^ sites are 
constrained in the following way: sites bl and b2 have the same a and b components but 
opposite c components. The same relation holds for sites b3 and b4. However, there are 
no relationships imposed by symmetry between the two set of sites mentioned above, bl/b2 
on one hand and b3/b4 on the other. For the Mn^+ sites, the moments are constrained to 
lie in the a6-plane but can have any magnitude and orientation within it. In the work of 
Munoz et al., the additional physical constraint of equal moments, not imposed by symmetry, 
was introduced. Here, the large number of independent magnetic reflections collected (204) 
allows the simultaneous refinement of all fourteen parameters. The result of the refinement 
is shown in Fig. [1] and the corresponding magnetic structure is displayed in Fig. [9] and 
Fig. [ini A R(F) factor of 4.5% has been obtained with a value of 18.2. BiMn205 has a 
very similar magnetic structure to that described previously for YMn205, for instance the 
antiferromagnetic ...-Mn^"'"-Mn^"'"-Mn^+-Mn^+-... chains in the a6-plane. This is at variance 
with the structure proposed by Munoz, in which there is a non-collinear arrangement of Mn^+ 
and Mn^"*" moments within one AFM chain. The angular deviation of the AFM chains from 
the a-axis (~8°) is less pronounced than in the YMn205 and HoMn205 compounds. Also 
unlike YMn205 and HoMn205, we found almost zero z-components for the Mn^"*" moments. 
The magnetic point group is m2m and the magnetic structures of the two possible domains 
are related by inversion symmetry, as previously described, resulting in flipping the central 
AFM chain in Fig. d 



IV. ORIGIN OF FERROELECTRICITY 



In this section, we discuss the possible mechanism that will promote a ferroelectric state 
in the magnetically ordered phase. First of all it is worth noting that irrespective of the 
microscopic mechanism involved, the point group symmetry m2m derived from the magnetic 
structures of YMn205, HoMn205, and BiMn205, is consistent with ferroelectric order only 
along the 6-axis. Secondly, it is important to discuss the possible microscopic mechanisms, in 
the light of the additional weak magnetic components observed along the c-axis for YMn205 
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and HoMn205 and not identified previously by powder diffraction. This small out of plane 
component introduces a modulation that resembles a cycloid, even though it does not corre- 
spond to a homogeneous rotation from site to site, as seen in Fig. |3]for example. However, 
this modulation could in principle give rise to ferroelectricity based on antisymmetric ex- 
change interaction, a mechanism proposed in several multiferroics, and its strength will be 
proportional to the spin-orbit coupling, potentially large for octahedral Mn^+. Obviously, 
the absence of such a modulation in BiMn205, multiferroic with a value of the electrical 
polarization comparable to YMn205 and HoMn205, does not support this picture. There- 
fore we argue that our model based on symmetric-exchange, and proposed initially in 
is the principal mechanism driving ferroelectric order. The recent observation of electro- 
magnonsj28] active only for certain polarization directions, also supports a model based on 
symmetric-exchange interaction. This model is relevant since the magnetic configuration in 
the a6-plane for the three compounds studied here is frustrated, i.e. the magnetic energy is 
invariant by flipping any AFM chain in the lattice. Such a configuration will induce atomic 
displacements on some of the ions in the crystal (we speculated that mainly Mn^"*" ions are 
involved based on structural anomalies related to this site) that will lower the magnetic 



energy 



291 . This comes at the cost of elastic energy, quadratic with the displacement. 



but for small displacements the gain in magnetic energy (linear with displacements) over- 
comes the elastic energy cost. Magneto-striction due to symmetric exchange should be much 



stronger than the DM interaction, which is a relativistic effect |l4j.ll5l.l30[|. However, it is pos- 
sible that this later mechanism contributes to the total electric polarization in the YMn205 
and HoMn205 compounds but should be of minor importance. In fact, the two mechanisms 
are likely to be strongly coupled. It is important to note that applying the inversion symme- 
try operation not only flips one of the AFM chains but also changes the rotation direction 
of the cycloidal modulation. Therefore the direction of the electric polarization in domain 
two is opposite to that in domain one, whatever the mechanism considered. The observation 
of cycloidal structures for many of these compounds is intriguing, and is unlikely to be a 
coincidence. In general, two mechanisms are known to generate helicoidal-type structures; 
competition between nearest and next-nearest neighbor interaction and direct DM interac- 
tion, resulting from the loss of a center of symmetry relating magnetic ions. It is easy to 
show that in none of the magnetic structures described here does the cycloidal modulation 
lower the next-near neighbor magnetic energy, regardless of the sign of the interaction. This 
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suggests that the observed c-axis modulation may be the effect, rather than the cause, of 
the loss of the center of symmetry relating Mn^"*" ions. 



V. SUMMARY 



The magnetic structures in the ferroelectric/commensurate magnetic regime of three com- 
pounds of the family i?Mn205 with R= Y, Ho, and Bi, have been determined by single 
crystal neutron diffraction. For YMn205 and HoMn205 the main features of the magnetic 
structures are in agreement with previous models derived from powder neutron diffraction 
characterised, in particular, by coUinear antiferromagnetic zig-zag chains in the a6-plane. An 
additional small component of moment parallel to the c-axis, and modulated in quadrature 
with the major components has been identified. This component was not detected in earlier 
experiments. The data also allow a more precise determination of the directions of the Ho 
magnetic moments. These results indicate the importance of super-exchange interactions in 
the magnetic ordering of Ho. The magnetic structure of BiMn205 is very similar to that of 
the Y and Ho analogs, but with the magnetic moments more closely confined to the a6-plane. 



2l|. 



It also contains zig-zag AFM chains, in disagreement with what was previously reported 
This comparative study strongly supports symmetric exchange as the principal mechanism 
leading to ferroelectricity since the small non-collinearity of the magnetic moments within 
chains along c is observed only for the YMn205 and HoMn205 compounds. This does not 
preclude a potential, but much weaker contribution from Dzyaloshinskii-Moriya interactions. 
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VII. NOTE ADDED IN PROOF 



During the final elaboration of this paper, we became aware of a recent manuscript 
by H. Kimura et al, which is now in print [31]. Kimura et al. investigated the crystal 
and magnetic structures of YMn205, HoMn205 and ErMn205 by single crystal neutron 
diffraction. Although the data analysis methodology is different from ours (they employ 
a large supercell rather than Fourier components), the agreement with our results on the 
compounds common to the two studies is excellent. Our work includes the refinement of the 
BiMn205, which, we believe, is crucial in clarifying the correlations between the different 
spin component and the electrical polarisation. 
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FIG. 1: Results of the refinement for the commensurate magnetic structures for three compounds- 
YMn20^, HoMn20z and BiMn20^. The experimental structure factors (times a scale factor) 
are plotted against the calculated ones. The agreement factors- magnetic R(F) factors and y^- are 
shown. See text for details of the magnetic models. 
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FIG. 2: (Color Online) Crystal structure of RMn20^ projected in the a6-plane (left panel) and 
ac-plane (right panel). Green, red and blue spheres correspond to Mn^"^, Mn^~^ and R^'^ ions, 
respectively. Mn-0 polyhedra are shown with the same color scheme. The thin black line repre- 
sent the crystallographic unit-cell. Correspondence between the labels and atomic positions are 
described in Table lllli 
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FIG. 3: (Color Online) Magnetic structure of YMn20^ projected in the ab plane. The structure 
is shown in two unit-cells, marked by thin black lines, along the a- and 6-axis. The green and red 
arrows represent magnetic moments on Mn^+ and Mn^^ sites respectively. Corresponding Mn-0 
polyhedra are shown with the same colors. Blue spheres represent Y ions. The different sites are 
labeled according to Table [Til 
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FIG. 4: (Color Online) Magnetic structure of YMn20^ projected in the ac-plane shown within 
four unit-cells along c and two unit-cells along a. The projection is shown separately for magnetic 
sites belonging to the first AFM chain (left) and the second chain (right) (See text for details). 
The green and red arrows represent magnetic moments on Mn^^ and Mn'^^ sites respectively. 
Corresponding Mn-0 polyhedra are shown with the same colors. Blue spheres represent Y ions. 
The different sites are labeled according to Table |TI1 
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FIG. 5: Magnetic structure of YMn205 projected in the 6c-plane and showing only the Mn^"*" 
moments. The figure shows the smah hehcoidal modulation generated by the out-of-phase c- 
component. To highlight this weak modulation, the moment have been scaled by a factor five with 
respect to what represented in the other Figures. 
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FIG. 6: (Color Online) Magnetic structure of YMn205 projected in the a6-plane (left panel) and 
the ac-plane (right panel) for the second domain, obtained by applying inversion symmetry on the 
magnetic structure of the first domain (Fig. [31 H]) 
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FIG. 7: (Color Online) Magnetic structure of HoMn20^ projected in the ah plane. The structure 
is shown in two unit-cells, marked by thin black lines, along the a- and 6-axis. The green, red and 
blue arrows (left panel) represent magnetic moments on Mn^"*", Mv?^ and Ho'^'^ sites respectively. 
Corresponding Mn-0 polyhedra are shown with the same colors. For clarity, in the right panel only 
Mn^+ and Ho^^ (belonging to an ordered Ho layer) magnetic moments are shown respectively 
in green and blue color. Lines connecting Mn^'^ and Ho"^'^ represent interactions (See text for 
details). The different sites are labeled according to Table Hill 
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FIG. 8: (Color Online) Magnetic structure of HoMn20^ projected in the ac plane shown within 
four unit-cells along c and two unit-cells along a. The projection is shown separately for magnetic 
sites belonging to the first AFM chain (left) and the second chain (right) (See text for details). 
The green, red and blue arrows represent magnetic moments on Mn^"'', Mn^^ and Ho^^ sites 
respectively. Corresponding Mn-0 polyhedra are shown with the same colors. The different sites 
are labeled according to Table IIIB 
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FIG. 9: (Color Online) Magnetic structure of BiMn20^ projected in the ah plane. The structure 
is shown in two unit-cells, marked by thin black lines, along the a- and 6-axis. The green and red 
arrows represent magnetic moments on Mn^+ and Mn^+ sites respectively. Corresponding Mn-0 
polyhedra are shown with the same colors. Blue sphere represent Bi ions. The different sites are 
labeled according to Table llVi 
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FIG. 10: (Color Online) Magnetic structure of BiMn20^ projected in the ac plane shown within 
four unit-cells along c and two unit-cells along a. The projection is shown separately for magnetic 
sites belonging to the first AFM chain (left) and the second chain (right) (See text for details). 
The green and red arrows represent magnetic moments on Mn^'^ and Mv?^ sites respectively. 
Corresponding Mn-0 polyhedra are shown with the same colors. Blue sphere represent Bi ions. 
The different sites are labeled according to Table IIVI 
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Symmetry operation 


Domain generated 


{1|000} 


Dl 


{T|000} * 


D2 


{m^2/o|000} * 


Dl 


{2.1000} 


D2 


{moy^lO^O} 


Dl 




D2 


{2,.ooliOO} * 


D2 


{2o^o|0i0} * 


Dl 



TABLE I: Action of the symmetry operations of the paramagnetic group on the magnetic structure 
of domain 1. Domain one and two are abbreviated by Dl and D2. Symmetry operations are shown 
in the Seitz notation. The * symbol indicates combination of the symmetry operation with complex 
conjugation operation. 
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a2 
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0.55(9) 


-0.50(12) 


0.125 


a3 


Mn^+ 


(0.58805 U.64921 U.5) 


3.01(5) 


-0.67(9) 


0.56(13) 


0.125 


a4 




(U.91195 U.14921 U.5) 


2.99(5) 


0.65(9) 


-0.60(13) 


0.125 



TABLE II: Magnetic Fourier coefficients obtained from least-square refinements of the single crystal 
diffraction data for YAIn205 at T = 25 K (Sec text for details). The Fourier coefficients for each 
site in the primitive unit-cell are given along crystallographic direction (x,y,z) and labeled M and I 
when real and imaginary, respectively. The phase for each site (Ph.) is given in units of 27r. Error 
bars are given for all parameters within parenthesis. 
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TABLE III: Magnetic Fourier coefficients obtained from least-square refinements of the single 
crystal diffraction data for B.oMn20^ at T = 25 (See text for details). The Fourier coefficients 
for each site in the primitive unit-cell are given along crystallographic direction (x,y,z) and labeled 
M and I when real and imaginary, respectively. The phase for each site (Ph.) is given in units of 
27r. Error bars are given for all parameters within parenthesis. 
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TABLE IV: Magnetic Fourier coefficients obtained from least-square refinements of the single 
crystal diffraction data for BiAln20rj at T = 10 X (Sec text for details). The Fourier coefficients 
for each site (M, real) in the primitive unit-cell are given along crystallographic direction (x,y,z). 
The phase for each site (Ph.) is given in units of 2tt. Error bars are given for all parameters within 
parenthesis. 
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